On Berezinskii-Kosterlitz-Thouless Phase Transition in 
Quasi-One Dimensional Bose-Einstein Condensate 

Vivek M. Vyas,^ Sandeep Gautam,^ and Prasanta K. Panigrahi^ 

^Indian Institute of Science Education and Research 
(USER)- Kolkata, Mohanpur, Nadia-741252, Indicia 
"^Physical Research Laboratory, Navrangpura, Ahmedahad-3 80009, Indi^ 

(Dated: January 26, 2013) 

Abstract 

We show that quasi-one dimensional Bose-Einstein condensate under suitable conditions can 
exhibit a Berezinskii-Kosterlitz-Thouless phase transition. The role played by quantized vortices 
in the two dimensional case, is played in this case by dark solitons. We find that the critical tem- 
perature for this transition lies in nano Kelvin range and below, for a wide range of experimentally 
accessible parameters. It is seen that the high temperature (disordered) phase differs from the low 
temperature (ordered) phase in terms of phase coherence, which can be used as an experimental 
signature for observing this transition. 
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The continuum version of 2D XY model appears in diverse areas of physics ranging from 
spin systems to superfluids jl-4|, with its Lagrangian being given by £ = — p(V6')^. It is well 
mown that this model exhibits singular vortex solutions: V0 oc ^ (in polar coordinates) 
l|, I2I]. Identifying two dimensional Cartesian plane with complex plane, allows one to write 
vortex solution in a rather elegant form : 

9{x, y) = 9{z = X + iy) = ±Imln(z — zq), (1) 

where Zq is the location of vortex. It is evident that the above solution exhibits a discon- 
tinuous jump as one traverses along a closed loop enclosing z^. This multivaluedness or 
discontinuity is a characteristic feature of these vortices and distinguishes them from other 
solutions. It has been shown that the Helmholtz free energy of a system of K vortices is 
given by 

F = K {2TTp-2kT)\n(^] , (2) 



where k stands for Boltzmann constant, L represents system size, and b is vortex core radius. 
Hence, one finds that below critical temperature Tbkt = free energy is minimized by 
having K = 0, which means that at temperature low enough vortices are thermodynamically 
unstable. However, when T > Tbkt, free energy minimization requires having K as large as 
possible. So the critical temperature Tbkt marks the occurrence of spontaneous proliferation 
of vortices, which makes condensate inhomogeneous. This is the celebrated Berezinskii- 
Kosterlitz-Thouless (BKT) phase transition 6|, |7|. This phase transition was found to be of 
infinite order, and unlike second order phase transition, here the ordered (low temperature) 
phase differs from disordered (high temperature) phase not in terms of symmetry of ground 
state but rather in terms of topology of the ground state ll,l2l]. In the low temperature phase, 
system exhibits a quasi long range order, where (e^^^^^e*^'-^^) ~ \x — y]'"' for some constant 
a. Above critical temperature the system has vortex excitations, which cause rapid phase 
variation, and one sees that correlation function ^e^^'^^^e*^^^-') ~ e~^'^~^' (/3 is a constant). 
This means that the quasi long range order present at low temperature is destroyed due to 
vortices, resulting in rapid decay in correlation. 

It is established by now that in the mean field approximation, physics of Bose-Einstein 
condensate (BEC) is well described by the Gross-Pitaevskii (GP) equation Sj: 

(-|^V^ + Vir) + f/ol^r) t) = (3) 



Here V{r) is the external trapping potential, and Uq = Airfi^a/m is the coupling constant 
for inter-atomic interaction, where m and a are atomic mass and s-wave scattering length, 
respectively. The above partial differential equation is an equation of motion describing 
space-time dependence of BEC, which itself is a ground state of a system of interacting 
bosons. In case when trapping potential is tightly confining along Z-direction and absent 
(or very weak) along the other two spatial dimensions, dynamics of BEC occurs only along 
the latter two spatial dimensions. The condensate function can then be written as \I^(r, t) = 
f{z)(f){x, y, t), where f{z) is a localized Gaussian wave packet along Z-direction, and y, t) 
describes condensate profile in the XY plane. Equation governing dynamics of such quasi-2D 
BEC can be obtained by integrating out degrees of freedom along Z-direction. Condensate 
density, in absence of any potential acting along XY directions, can be expected to be 
stationary (temporally) and spatially uniform: \(j){x,y,t)\'^ = p, where p = constant, and </> 
field can be written in terms of the polar variables: = yfpe^^^^'^\ It can be shown that in 
static case, Lagrangian for quasi-2D BEC reads as oc —p(V9Y, which exactly matches 
with that of 2D XY model. Hence, one finds that dynamics of quasi-2D BEC is actually 
captured by XY model, and so one expects to observe BKT transition in this system. 

Above discussion is oversimplified, and in reality, bosons are trapped not only in Z- 
direction but also in XY directions. Because of this, there appears a confining potential 
defined over XY plane. This makes the system inhomogeneous, and hence translational 
invariance is lost. In the last few years, there have been significant activity, on both theoret- 
ical and experimental fronts, to gain more understanding about nature of BKT transition in 



these systems |9Hlli. A 



ew years back, Hadzibabic et al, reported a direct observation of 



12| . Further, in recent years study of dynamics of vortices in BEC 



BKT transition in BEC 
systems especially, vortex dipoles (vortex-antivortex pairs) has received significant attention 
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18| . Vortex dipoles are also produced as the decay products of dark solitons in 2D BECs 



due to onset of the snake instability |l9 |. 



In light of above discussion, it is natural to ask, if a BKT phase transition occurs in BEC 
in spatial dimensions other than two. The answer is crucially dependent on another question: 
Do the counterparts of vortices exist in dimensions other than two ? Analogues of vortices 



in three dimension are well known, they are vortex tubes or rings 
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21| . Being spatially 



extended o 
nontrivial 



ejects, interaction between them is complicated, which makes the problem rather 
23| . Instead, we ask whether there are any analogues of vortices in one 
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dimension? If yes, then does one dimensional BEC display BKT type phase transition? 
Answer to the first question is known. The so called dark soliton solutions are known to 



be analogues of vortices in one dimension 2J, |25[. In what follows, we answer the second 



(which is also used in the 



question affirmatively. After introducing dark soliton solution of GP equation, we indicate its 
analogy with two dimensional vortex. Thereafter, we show using exactly similar free energy 
argument as was originally used by Kosterlitz and Thouless 
above discussion) that it is indeed possible to have a finite temperature BKT phase transition 
in quasi-one dimensional BEC. We estimate the critical temperature for the same and find 
that it can be tuned over a broad range by changing scattering length and number density. 
Subsequently, we also show that the long range phase coherence exhibited by the system at 
low temperature vanishes at the temperatures above critical temperature. We conclude by 
mentioning that it may be possible to experimentally realize this phase transition. 

As mentioned earlier, GP equation (|3]) governs motion of order parameter \E'(r, t), which 
satisfies normalization condition J |\E'(r, )f:)p(ir = A^, where A^ is the number of the bosons. 
In present case, we consider trapping potential to be harmonic along the radial direction 
and infinite square well along the axial direction, i.e., 

{m r, 2 -r L L 
— UJ r if < z < 
2 \ ^ (4) 
oo if \z\ > — , 
I I - 2' 

where u is the radial trap frequency and L is the size of infinite square well potential along 
axial direction. The infinite square well potential can be created by applying two blue 
detuned laser beams with narrow beam waist and sufficiently high intensity at \z\ = L/2. 

Since the potential becomes infinitely high at and after \z\ = L/2, condensate density 
vanishes at \z\ = L/2 and increases smoothly as one moves towards the center. The char- 
acteristic length scale over which BEC recovers its bulk density tiq is called healing length 
or coherence length f = The characteristic length of the BEC along radial direction 

is oscillator length Oosc = a/^/ (mu). When a^sc ^ ^, Qosc ^ L, and BEC is in the weakly 
interacting regime 2'iraN J r|\E'(r, z)p(ir <^ 1, the order parameter can be factorized into 



radial and axial parts 26 1 



^ir,z,t) = xir)^iz,t), (5) 

where x(^) is the normalized ground state of radial trapping potential, which can be con- 
sidered to be time independent provided, kT -C hw. We consider x(r) to be normalized to 



unity {2tt J r\x{r)\'^dr = 1); then from Eq.Q and Eq.(l5]), after integrating out the radial 
order parameter, GP equation reduced to one dimension reads, 



(6) 



2m dz"^ 

where u = 2ahw and V^lz) = hw for —L/2 < z < L/2 and zero otherwise. The order 

L/ /2 

parameter satisfies normahzation condition: / i ,2dz\ilj(z)\'^ = N. It is well known, since 
the pioneering work of Zakharov and Shabat [271], that equation ([6]) allows for grey solitons 
solutions. In terms of atoms per unit length along axial direction a{z) = \iIj{z)\'^ and a 
dynamical phase, solution at t = is given by 



tan 



ao 1- 



2$ij{z) = y/a{z)e'^^^\ where 

COs2(^) 



cosh {zcos{9)/Q 



and 



■ a/(1 — u^/cl/) tanh 



v/(i-uV4)/ 



(7) 



c 

Here velocity of soliton is given by u and velocity of sound is c^. Parameter 9 is defined as 
6 = sin~"'^(M/c^), whereas ( is defined as C = '^^{^q), and ctq = UQiial^^ is the linear density 
far away from the soliton, i.e., in the domain v oc (f)'{z) 0. Notice that above soliton 
becomes 'dark', i.e, condensate density a{z) vanishes, at 2; = only when u = 0. Therefore, 
by construction, a dark soliton in this model is static. Amazingly, in the limit u approaches 
zero, from Eq.Q, one sees that the phase of the dark soliton behaves like a Heaviside step 
function with height vr and becomes discontinuous at z = ,i.e., 

TT 





TT 

2 



if 2 < 0, 

if z = 0, 
if 2 > 0. 



It is this discontinuity that makes the dark soliton a ID counterpart of vortex. A cautious 
reader may ask that the above mentioned solitons are strictly solutions of Eq. (|6]), when 

30X tj^e 



2J], that 



^{z) is defined over whole of real line. In present case, because of particle-in-a- 
potential, it is only defined in domain ^ < z < ^. It was shown by Carr et. al. 
dark soliton is indeed a genuine solution in this case too. It can be easily shown that the 
en ergy required to create a dark soliton in otherwise homogeneous and uniform condensate 
is 



4^2 



^(o^aoaosc^- 



In terms of (7(2), the criterion for applicability of Eq.([6]), i.e., 27ra J r|\l/(r, z)p(ir ^ 1, 
becomes o'{z)a ^ 1. The transition between weakly interacting to strongly interacting high 
density limit occurs for a{z)a ~ 0.25 [29]. If is large (A^ — > 00) such that ^/L ^ 0, i t 



can be argued that a dark soliton can be created without disturbing the other soliton 



24|. 



^n reality, however, is finite, and ^/L quite small but non-zero. It was shown in Refs. 
and we have also checked numerically that above argument still holds as long as the 



distance between two dark solitons is much greater than the healing length. So the energy 
required to create K solitons is Ek = Ke. Entropy associated with K solitons is given by 
KkhiW , where W is number of independent states soliton can occupy, and it is easy to see 
that W = L/^. Hence, Helmholtz free energy for K solitons is: 

F = K ^^y^ V (o-oa)o-oaosc^ - kT In , (8) 
which changes sign at a characteristic temperature 

Tbkt = 3fc hi'^L/O (^o'^)^o'^°«^^- (^) 

So exactly like in 2D BEC, we see that in ID BEC also there exists a finite critical temper- 
ature Tbkt, above which a uniform homogeneous condensate becomes thermodynamically 
unstable and spontaneous proliferation of dark solitons takes place. This indicates that ID 
BEC exhibits a non-symmetry breaking BKT phase transition. Note that critical temper- 
ature obtained above, depends on various tunable parameters like a and ctq- In a typical 
cold atom experiment, the maximum density in trap (no) can vary from 10^^ — lO^^cm"^ 
30| . whereas typically u varies from a few Hertz to kilo Hertz. For u = 130 Hz, L = 500^, 
density around lO^^cm^^, and scattering length around 50 oq {ao is Bohr radius), one finds 
that Tbkt lies in sub nano Kelvin range (~ 0.1 K), which is an order of magnitude less than 
a typical BEC transition temperature. Figure ([1]) shows dependence of Tbkt as a function 
of scattering length and density, and it is clear that it can be conveniently tuned by varying 
density and/or scattering length at will. 

Dimensionally reduced GP equation can be obtained by extremizing the action S = 
J dxdt^, where 



^ dt 2m 



dip ^ 



dz 
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In case when background condensate density is constant, only dynamical degree of freedom 
left is phase, since ip{x,t) = ^yooe*^*^^'*-*. Dynamics of 9 field is governed by Lagrangian 
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FIG. 1. Density plot depicting Tbkt (in nano Kelvins) as a function of bulk density no and s-wave 
scattering length a 



Provided that there are no dark solitons, static spatial phase correlation function can be 
easily found from the above Lagrangian and is 

(^(^)^(0)) = 

which clearly indicates that system indeed has phase coherence, a characteristic feature 
of any BEC. However, situation changes when dark solitons are present. Phase field now 
can be decomposed into a regular continuous part and a discontinuous part, i.e., 9{z,t) = 
6reg{z,t) + 9ctis{z), respectively. Substituting this in the above Lagrangian, one finds 



dt ' 2m "''^^ dz'^ 

Teg 



+ — ^diS r. o 1 VdiS 



2m oz-^ m oz'^ 

The effect of dark solitons on phase correlation can be inferred by first integrating out the 
singular part of field to yield an effective action for regular field. Very interestingly, one finds 
that integration of Odis gives rise to a term — ^f^^reg^^?^, which exactly cancels the gradient 
term already present in the action. Hence, phase field in presence of solitons does not 
propagate, and system losses long range phase correlation. In the absence of solitons, phase 



variation along the condensate profile is rather slow, and hence such a condensate would 
yield a well defined interference pattern, an authentic signature of coherence. However, 
when solitons are present, phase receives random discontinuous kicks along the profile, as 
a result phase ultimately gets averaged out. Such a condensate would certainly exhibit 
poor interference pattern, an indication of loss of coherence. Hence we conclude that high 
temperature solitonic phase is actually a disordered phase, and transition from uniform 
condensate to solitonic state is a phase transition. 

We have shown that like 2D BEC, ID BEC can also exhibit a BKT phase transition at 
finite temperature. The role played by vortices in 2D is, in ID case, played by dark solitons. 
It is found that critical temperature for this transition depends on atom number density 
and scattering length, both of which can be finely controlled in a typical experiment. It 
is this particular feature that could facilitate and strengthen possibility of an experimental 
realization of this transition. One may wonder whether above discussed phase transition 
indicates violation of Coleman-Mermin-Wagner-Hohenberg theorem This theorem 

states that a continuous symmetry in given infinite system (not bounded by a finite vol- 
ume), which is homogeneous and isotropic with short range local interactions, can not be 
broken at zero temperature in one spatial dimension and at finite temperature in two spa- 
tial dimensions. The proof of above theorem is based on the fact that in such a system, 
breaking of a continuous symmetry results in existence of gapless Goldstone modes which 
by construction lead to infrared divergence. It was shown that in one spatial dimension, this 
divergence is so severe that it does not allow ordered phase to exist, and hence a continuous 
symmetry is not broken even at absolute zero temperature. As is obvious, this theorem 
does not hold for our system of ID BEC for several reasons. Firstly, we are dealing with a 
system which is externally confined into a finite volume by trapping potentials, and hence 
is neither homogeneous nor isotropic and certainly not of infinite extent. Further, we are 
actually looking at a quasi ID BEC, where BEC itself lives in 3D but at sufficiently low 
energies, degrees of freedom along other dimensions are frozen. Moreover in BKT phase 







transition, no symmetry of Hamiltonian is spontaneously broken by ground state 
So we see that Coleman-Mermin-Wagner-Hohenberg theorem does not hold for our system. 
Possibilit y o f BKT transition in ID systems (in thermodynamic limit) have been studied 



in Refs. 
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32] . and both find transition temperature to be absolute zero. Interestingly, 



taking L — )• oo in equation ([9]), one sees that Tbkt goes to zero, which agrees with their 



results. Issues pertaining to critical behavior and universality class of this transition are 
currently being investigated and will be published in due course. 
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